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ABSTRACT In this paper, a co-infection of malaria and visceral helminthiasis with immune 
stimulation and impairment is studied. We have assumed that an individual gets malaria in- 
fection during invasion by helminths larvae. In absence of immune response, our results show 
that the antigens invade the blood system if the rate of red cell rupture per invading merozoite 
is greater than one. If fewer merozoites are released, the initial state is globally asymptotically 
stable. If more merozoites are released, there exists malaria-only endemic point. However, 
both antigens coexist if the mean infection burden is greater than one. In this case, there is 
a threshold value for drug action below which no recovery of host is expected. In presence 
of immune response, three equilibrium states exist. The initial invasion state, and secondly, 
the unstable state when the immune population has been triggered and the antigens have 
been eliminated. The third state represents the endemic state which is stable if the infection- 
induced rupture rate per erythrocyte with immune activation is less than the total mortality 
of the erythrocytes. A model for severity of the co-infection shows that an immune response 
will be delayed until immunological barrier values for malaria and helminths are exceeded. 

Keywords: Within-host co-infection; Malaria; Visceral helminthiasis; immune 
response, suppression, impairment; permanence; persistence. 



1 Introduction 

Specifically, in this model of immune reaction, we explore if an asymptotic decrease of antigen quantity is 
reached, the cause of periodic course of the illness, and if there is unlimited growth of the antigen quantity 
within the human, and under what conditions. The model allows a description of threshold relationship 
between the infection process and the initial dose of the antigen. We seek to determine whether destruc- 
tion of either malaria parasites or helminths larvae is possible without spread of infection, destruction 
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after either infection, coexistence with specific antibodies, a possibility of recurrent course of infection of 
either disease, and if there is unlimited multiplication of either malaria parasites, or helminths larvae. It 
is expected that this co-infection model will give a new insight to researchers in designing experiments for 
discovering the underlying mechanism of immune response and suppression in within-host malaria and 
visceral helminthiasis co-infection. We have assumed a normal functioning of the immune system where 
there is no distinction between the cellular and humoral components of immunity fighting the malaria 
parasites and helminths larvae that have penetrated into the human blood circulation but that a human 
being has such defense components, we refer to as antibodies no matter whether it is a cellular-lymphoid 
system, or humoral- immunoglobulin system of immunity [11]. 



2 Model formulation 

The biologically known within-host dynamics of visceral helminths are such that when a mosquito or 
fly carrying microfilarial larvae bites a new host, it deposits larvae into the subcutaneous tissue of the 
host. The larvae find their way into the lymph vessels where they take about a year to mature. The 
mature adult female produces larvae called microfilaria. These escape from the lymphatics to reach blood 
circulation. It is from here where mosquitoes biting the infected person ingest them [11]. The ingested 
microfilaria develop in the mosquitoes to reach the infective stage of microfilaria and then transmit them 
to the next human victim when they bite him. This is similar to the transmission of malaria parasites. 
The injected larvae migrate from blood stream into the lymphatics thus completing the life cycle; when 
they mature into adult filarial worms. Microfilaria appear in plenty in blood circulation at distinct times 
of the day, usually in the evenings and at night. This characteristic phenomenon is referred to as period- 
icity and appears to be related to the biting habits of the transmitting mosquitoes. 

In humans, malaria parasites grow and multiply first in the liver cells and then in the red cells of the 
blood. In the blood, successive broods of parasites grow inside the red cells and destroy them, releasing 
daughter parasites (merozoites) that continue the cycle by invading other red cells. 

The blood stage parasites are those that cause the symptoms of malaria. When certain forms of blood 
stage parasites (gametocytes) are picked up by a female Anopheles mosquito during a blood meal, they 
start another different cycle of growth and multiplication in the mosquito. After 10-18 days, [11], the 
parasites are found (as sporozoites) in the mosquito's salivary glands. When the Anopheles mosquito 
takes a blood meal on another human, the sporozoites are injected with the mosquito's saliva and start 
another human infection when they parasitize the liver cells. Thus, the mosquito carries the disease from 
one human to another (acting as a vector). Differently from the human host, the mosquito vector does 
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not suffer from the presence of the parasites [11]. 

Consider a population of N e (t) erythrocytes in an individual. Let the uninfected erythrocytes be U e (t). 
Assuming an invasion of malaria parasites (merozoites) with concentration M(t). At the same time, as- 
sume that the individual is attacked by a mass of helminths larvae H(t). Each of these two parasites 
invade random red blood cells and imprcgnants them. When the invasion is successful, the uninfected 
cell becomes infected. Let the body of infected red blood cells be I e (t). Therefore N e (t) = U e (t) + I e (t). 
After infection, immune responses against the malaria parasites and helminths larvae are triggered. A 
wide body of evidence indicates that protective anti-blood stage immunity is dependent upon CDA + T 
cells [44]. Therefore, effective immunity is considered as a population of activated T cells and assume 
that the resulting immune response is a direct function of their density. The magnitude of the immune 
response is proportional to the density of the immune cells. The immune cells augment the clearance of 
the merozoites, helminths larva, and infected red blood cells from the body. Wc assume that immune 
activation is proportional to the density of the infective stages and the precursors are not limiting. Define 
T(t) as the number of antibodies recruited when their resting precursors contact free merozoites, infected 
cells and free larvae at net rates 7 M M, 7 f I e , and j H H respectively. Let the rate at which the immune 
cells expand be p. This rate of expansion encapsulates the positive feedback upon the immune system. 
We further assume that a regulatory negative feedback force operates to suppress immune population 
growth at a rate proportional to the square of its density, bT 2 [22]. This function implies a regulation of 
the response at high antigen concentrations, since convex or plateauing relationships between the rate of 
T cell proliferation and antigen concentration is observed [24, 30]. Thus, we have clearance rates of free 
merozoites and larvae due to B cells and macrophages given by XmMT and XrHT respectively. The 
clearance of the infected red blood cells due to T cells is given by XjI e T. The rate of change of density of 
immune cells is described by their proliferation and deaths rates. They proliferate in response to contact 
with free-merozoites, free larvae and infected red blood cells at rates 7 M M, "f H H and jjl e respectively. 
We define the natural clearance rate of the antibodies as /i t . All parameters of the model are positive 
real numbers. 

Let the uninfected erythrocytes be recruited from the bone marrow at a per capita rate A. We further 
define removal of the erythrocytes due to aging as a per capita rate \i e . Let olm be the per capita rate of 
erythrocyte invasion and comprises of the dual probability of contact between an uninfected erythrocyte 
and a merozoite, and of such a contact resulting in productive invasion. The corresponding invasion term 
for helminths larvae is an- The infected red blood cells experience a per capita death rate S M due to 
infection- induced rupture where 6 M is large compared to fi e . Each ruptured infected cell produces an 
average of r merozoites such that the net rate of merozoite production is rS M I e which is reduced by 
drug action effect e. Free merozoites die at natural rate fi M . If the helminths larvae are not ingested 
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by a feeding fly or mosquito, they die at a per capita rate /j, h . The larvae suffer additional death due 
to prophylaxis at an average rate S H . If we assume that merozoites and larvae invade erythrocytes at 
random, that is, have no predilection for erythrocytes of a particular age or type [22], then the above 
definitions and assumptions lead to the following set of non-linear ordinary differential equations: 



^ = A - n e U e - a M U e - a H U e 



= a M U e +a H U e - n e I e - 6 M I e - X T I e T 



rS M (1 - e)I e - m m M - a M U e -\ M MT (1) 



dt 

dM 
dt 

<f= U-^ H H-6 H H-a H U e -X H HT 

f = lM M + ll I e + lH H+pT-bT 2 -^ t T 

First we analyze this system for the case when there is no immune response and determine the 
existence and nature of the stationary states. 

2.1 Criteria for invasion and persistence in the system without immune re- 
sponse 

As in [28], without immune response System (1) simplifies to 



dU e 
dt 



A — jj, e U e — a M U e M — a H U e H 



= a M U e M + a H U e H - jj, e I e - S M I e 

(2) 



dt 

dM 
dt 

<f= n- n H H -S H H -a H U e H 



rS M (1 - e)I e - ji M M - a M U e M 



Using the next generation matrix method we show that the basic reproduction number is given by 



—a M rS M {l - e) 

Ro = it ( 3 ) 

(fie + S m )(hm + —a M ) 



Next we determine the nature of stability of the system. The first steady state is when there is no 
invasion by malaria parasites given by U* = — , I* = M* = 0, H* =11, II <C 0. The endemic 
state is too complicated to obtain useful analytical results. We use numerical simulation to integrate 
numerically and obtain the solutions as shown in Figure (1). These graphs give an example of how 
endemic stability is attained. The phase portrait indicates that the system moves towards a fixed point 
at endemic equilibrium. 
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Figure 1: Graphs for the evolution of uninfected and infected erythrocytes with time (A). In the graph, 
the circles represent uninfected cells while the crossed line is for infected cells. The phase space portrait 
in the uninfected U e and infected I e space is shown in (B), while the evolution of malaria M and helminth 
larvae H with time is shown in (C), here, the diamond line is for malaria parasites while the star line 
represents the helminths larva, e = 0.5, A = 250, ^t e 0.025, (ih = 20, [Im = 48, 8m = 1,6h = 1,&m = 
0.025, a H = 0.005, r = 16, II = 25, -y M = 0.1, 7/ = 0.1, j H = 0.1, p = 0.1, b = 0.000000001, fi t = 0.05, A/ = 
0.00000001, A M = 0.00000001, \ H = 0.00000001. 



Lemma 2.1.1 In the absence of immune response the malaria parasites and helminths larvae are only 
capable of initially invading the blood system if criteria equation R D > 1 . Biologically, if the rate of red 
cell rupture per invading merozoite is greater than one. 
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Next we determine the local stability of the equilibrium states. Without immune response, the Jacobian 
of the system at endemic equilibrium is given by 



J{E") 



"Me ~ <*M M * ~ a H H * 

a M M*+a H H* 
-a M M* 





~ S M ~ Me 

-a M M* +r5 M (l-e) 







The malaria invasion-free steady state is stable if 



(4) 



R 2 = a M rS Ml^( 1 - ( -) <1 

+ S H + a H lt^ S M +Me+M M +a M + + Me)(M M +" M £)j 



(<5 M +Pc)(p M +a M ^) 



< 1 



(5) 
(6) 



Proof: We show how the local stability of the initial invasion state without immune response is obtained. 
The Jacobian of system (2) at endemic equilibrium is given by 



J(E* 



-/i e -a M M* -a H H* 
-a M M* 





S M - Me 
-a M M* +r8 M (l -e) 



When there is no invasion by either pathogen, U° 
invasion-free state given by 



(7) 



We therfore obtain the Jacobian at the 



-Me 







~ & M ~ Me 



r5 M (l-e) -,, 




1 Me 
A 



M M , 



n Me 

_A_ 

Me 



H H "ff 



(8) 



with distinct eigenvalue — fi e . The remaining three are obtained from the 3x3 square matrix given by 



~K - Me 
rS M (l-e) 




M Me 



"Mm - a, 







-M H ~S h -a H — 



(9) 



with characteristic polynomial 



A 3 + a x \ 2 + a 2 X + a 3 = 



(10) 



6 



where 



ffli = 5 M + n e + n M + m h +S h + (a M + a H )±) > 

«2 = (Mh + * H + "h + Me + Mm + »M £) + + Me) (Mm + »M^) [1 - . 

where i? 2 = "m^mAC 1 -') 

> since i? 2 < 1 

a 3 = (m* + K + a H ±){5 M + [i e )(p M + a M ±)[l - R 2 J, 

> if R 2 < 1 



(11) 



where R 2 



( S M+^)(M M + a M-t) 



ttt a i a 2 — ^3 > implies that 



>{l-R 2 ). (13) 



This is greater than zero (> 0) if and only if R 2 < 1 and R 2 < 1, where -0 = (<^ M + Me + M 



M ' M , 



A. 
/'< 

result is obtained: 



Mff + °~h + a ff ^ )• Therefore, all eigenvalues of the system have negative real parts. Hence, the following 



Lemma 2.1.2 The invasion-free state is locally asymptotically stable as long as R 2 < 1 and R 2 < 1. 

For global stability of the initial invasion state, consider the following Lyapunov function, 

C = I e + M + H 

which is positive definite V I e , M, H as t — > oo, with orbital derivative given by 



C'= I' e + M' + H' 

= a M U e M + a H U e H - /j, e I e -S M I e 

+rS M (l - e)I e - fi M M - a M U e M + 11 - fi H H - 5 H H - a H U e H 
-Hel e - 5 M I e + rS M (1 - e)I e - n M MM + II - n H H -S H H 

< ( rS M (1 - £ ) - Me - S M )I e - n M M - n H H -6 H H 
(Me + SjC^frP 1 - iKe - M M M -h h H- 5 h H 

( Me + 5 M ) (i? 2 - l)I e - Mm M - (p a +5 H )H, where i? 2 = 

< if < 1 

It is important to note that helminths larvae invade the blood stream a year after an individual is bitten 
by the fly. It is therefore biologically meaningful to assume that H — > as t — >• oo. With this assumption, 
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our system has a maximum invariant set for £'=0 if and only if Rq < 1 holds and there is no new invasion 
of the red blood cells. Therefore, all the trajectories starting in the feasible region where the solutions 
have biological meaning approach the positively invariant subset of the set where C— 0, which is the set 
where I e = M = and H -> 0. In this set N t 



-7- as t — > +00. This shows that all solutions approach 
the initial invasion steady state. Thus, when R% < 1, malaria and helminthiasis will be eliminated from 
the blood circulation. If R% > 1, then £> for I e , M, N e , H close to (0,0, ^,0) in T except I e = 0. 
Thus, the following conclusion is made: 

Lemma 2.1.3 The initial invasion state is globally asymptotically stable in T if R^ < 1. In this case, 
the co-infection of malaria and helminthiasis will be kept under control. This would mean that fewer 
merozoites will be released per rupturing infected erythrocyte. 

When Rl > 1 then V [r |l~ e ? ] > 1 implying that (1-e) > tti+hd. Thus, e < 1- { ±^m1. This gives a 

\t J -"+°M' T °M r °M 

percentage value of drug action that is sufficient for treatment of the co-infection of malaria and visceral 
helminthiasis in absence of immune response. Thus we have the following result: 

Lemma 2.1.4 e — 1 — ^^m1 j s ^ e threshold value of the drug action necessary and sufficient to treat 
the co-infection of malaria and visceral helminthiasis in absence of immune response. 

The biological feasible region for system (2) is the invariant simplex in the positive cone of i?* given by 
r = {(U e , I ei M, H) 6 R\_} including all of its lower dimensional boundaries. Mathematically, system (2) 
will be regarded as a system in R 4 ^ with an invariant manifold T of dimension 3. Re-writing (2) in the 
form 



x = f{x), z = ^{x)z, 



(14) 



'4^ 



V 3 / 

where z = (z\, Z2, Z3, Z4) G R = R . The third additive compound -i^- for system (2) using (4) is 
given by 

a/Pi 



dx 



■[Me + a M M + a H H + 5 M + fi e + n M + [a M + a H )U e ]I 



where is a matrix given by 



(15) 



a H U e 


a„H 



-a H H 





-H H H -5 H H + a M U e + fj,, 
-a M M + rS M (l-e) 



a M M 



-a H U e 



-a H U e 
a M U e 



5 m + Me- H H H -6 H H 
a M M + a H H p 
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(16) 



with p = —p. H H — 5 H H + /j, e + a M M + a H H. To show asymptotic stability of the endemic equilibrium, 
we use the method of first approximation [26]. Using the spectral properties of the second compound 
matrices, we state the following Lemma. 

Lemma 2.1.5 Let A be an n x n matrix with real entries. For A to be stable, it is necessary and 
sufficient that 

1. the third compound matrix A^ is stable, 

2. {-l) n det{A) > 0. 

Proposition 2.1.1 The endemic equilibrium is stable if the following inequalities are satisfied 



Ar6 M (l-e)I* e > 2(a M M*[/*) 2 
6 M W* e > 2(a H U:H*) 2 



(17) 



that is, if the mean infection burden is greater than 1. 

Proof: The third additive compound matrix of System (2) , , is given by (15) . For E* = (U* ,I*,M*,H*), 
and the diagonal matrix D = diag(£7* , I* , M* ,H*), the matrix J^(E*) is similar to DjW(E*)D~ 1 . This 
is given by 



-[Me + a M M* + a H H* + S M + p e + Mm + (a M + a H )U*]I + V 



(18) 



where ip is the following matrix 



a H u: 



a H H* U* 

M* 
-a H H*U c 
H* 





-H H H -6 H H + a M U: + 

(-a M M*+r6 M (l-<i))i: 
M* 



-a H U;M* 



H' 



(a M M*+a H H*)M* 



-a H U'*H* 



(19) 

with 9 = -n H H-6 H H + (j, e + a M M*+a H H*. The matrix J^(E*) is stable if and only if D {E*)D~ l 
is stable, for similarity preserves the eigenvalues. Since the diagonal elements of the matrix DJ^ (E*)D~ 1 
are negative, an easy argument using Gersgorin discs shows that it is stable if it is diagonally dominant 
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in rows, [45]. Set \i = {31,52,53,34} where 



9i = - W + a M M* + a H H* + S M + Me + Mm + a M U* + a H (M* + H*)} 

.92 = -[Me + a M M* + a H H* + S M + ^ + ttfl [/;] - Mh H - S H H + a ^ M * + ^mEL 

.93 = - [Me + a M M* + a H H* + fj, M + (a M + a H )U*] 

T M , T M , f^H 11 U H 11 

94 = - [S M + A*e + Mm + Om + a H ) C/ e1 



At steady state, System (2) can be written as 



|[^- Me + ( /Ue + ^ M )i] = (a M M* +<*„#*) 
JL-^ H H-5 H H =a H U; 



Therefore, substituting Equation (21) in (20) we obtain 



9i= - [f + 2% + Me(l + Jt) + <5 M [1 + 7 e *(^T + ^)] + a H (M* + H*) 

92 = - + + Me(l + 2%) + + 2^) - ^(M* + H") 

93= -[^ + ^ I + We + 5 M )§ ! + ^+a M I* e -'^E: 

94 = - [S M + 2^] + Me(l - + £ + a H U: + ^ e - ^] - 



This gives /x < which implies diagonal dominance as claimed and thus verifies the first condition. 
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From Jacobian 4 and Equation (21) , 



det[J(E*) 



-He -a M M* -a H H* 
a M M*+a H H* 





-$M - Me 

-a M M*+rS M (l-e) 



"Mm " a M U e 




~ a H U e 





-^llT* - Me + (Me +& M )jf*\ 





-<5 M - Me 

-a M M*+rS M (l-e) 
-a„H* 



~ a M U e 

M* 






n 



{ 



H 



-($M +Me) 



1 r A / c N J t 



rS M (l-e)i; 



M* 



)( 



M* 



(23) 



V u m ' Men 2H»U;M' ' 2H»M»U; ~ 2H*M'U* 
r6 M (l-e)^m: na^M-U: r5 M (l- £ )a 2 w u;i;H' . 



2H*M* 



AW> M (l-e)/ e * 



M* 



2a 2 M* 2 */* 2 

M c 

^ nj : _ 

2a 2 U? 2 H* 2 



- 1 



-} 

2H'M* 



(24) 



Note that in the hrst expression, ri5 " ^g^- 5 is the ratio of the net rate of mcrozoitc production r<5 M (1 — 
e)I e per ruptured infected cell in presence of drugs to the total number of free merozoites that invade 



susceptible erythrocytes a M U e M, while 



is the invasion rate per recruited erythrocyte. The second 



expression gives the mean infection burden of erythrocytes, while the third expression 2^~w^w^ 

e H e 

gives the net rate of helminths invasion of red blood cells per released larvae - — y, H , and the net rate of 
rupture of infected erythrocytes in the presence of helminths invasion - S fji° H , at endemic equilibrium. 
Therefore, det[J(E*)] > and this completes the proof. □ 



2.2 Criteria for invasion and persistence in the system with immune response 



Consider the case when there is immune response to the infections. In this case, we have three categories 
of equilibrium: the hrst is the initial invasion state when the number of helminths larvae released in 
the blood is small and hence the host is not displaying any symptoms and the immune response is not 
triggered. This is the naive equilibrium where U e = — and I e = M = T = 0, H « 0. Secondly when 
^ = 0, there is a second equilibrium point where U e = — and I e = M = 0, and H « but T = p ~^* . At 
this steady state, the immune system has been triggered and the malaria parasites and helminths larvae 
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have been eliminated and the system has returned to an infection-free state with a degree of residual 
immunity. The third type of equilibrium represents host-parasite co-existence and is characterized by 



non-zero levels of I ei M and H, U e < — and 

Me 



T > 



V-IH 



(25) 



The endemic state is too complicated to be useful. However, numerical integration shown in Figure (2) 
give an example of how this endemic state can be attained. We make the following conclusion: 

Lemma 2.2.1 In the "presence of immune response, there are three equilibria states. The first is the 
initial invasion state where the helminths larvae die off without being ingested by a feeding fly. The 
second steady state is when the immune population has been triggered and the malaria parasites and 
helminths larvae have been eliminated. In this case the system returns to an infection-free state with a 
degree of residual immunity. The third state represents non zero values ofU*, I*, M* and H* . 

From this result, we have: 

Lemma 2.2.2 In the absence of infection, the immune population (and therefore the naive equilibrium) 
is inherently unstable and will not return to zero once it has been triggered by the presence of malaria 
parasites and helminths larvae. 

From Figure 2, we note a delayed immune response and thus an increase in pathogen concentration. 
Once the immune response is triggered, there is a significant reduction in I e , M and H. 

Next we analyze stability of the states by introducing small perturbations around these states. The 
Jacobian of the System (1) is given by 



Je = 



■hi J\2 
•hi Jii 



(26) 



where 



-H e -ot M M* -a H H* 
a M M*+a H H* 



-S M - ^ - \iT* 



, Jl2 



~ a M U t 
a M U e 



-a„u: 





-\ M M* 



J21 



-a M M* -a M M* + rS M (l-e) 



li 



J22 





7m 



1h 



p-fit- 2bT* 
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Figure 2: Graphs for the evolution of uninfected and infected erythrocytes with time (A); the dashed 
line is for un infected erythrocytes, while the solid line represents the infected erythrocytes. In (B), the 
evolution of immune response with time is shown, while in (C), evolution of malaria M and helminth larvae 
H with time is given. In the graph, the dashed line malaria parasites while the solid line represents the 
helminths larva. Parameter values used are e = 0.5, A = 250, /i e 0.025, fit = 20, hm = 48, <5jvf = l>£ff = 
l,a M = 0.025, a H = 0.005, r = 16,11 = 25, 7 M = 0.1,7/ = 0.1, j H = 0.1, p = 0.1, b = 0.000000001, fi t = 
0.05, A/ = 0.00000001, X M = 0.00000001, X H = 0.00000001. 
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-(p - Ht) 



13 



(27) 



with distinct eigenvalues — fi e and —(p — f-it) since p>/i t . The remaining three are given by a 3x3 square 
matrix given by 



rMl-e) 




M Me 



"Mm - "a 



P-JM 



"Mi 



a, 



(28) 



with negative eigenvalues — \p, H + a H + Xh p and — ^ a + e + ^/a 2 + e 2 — 2ae + Idf where a = 
(<L + Me + A/ 2 ^), e = (/z M + a M A + Am^), = a M A and / = rS M (l - e) respectively. The last 



eigenvalue given by — | a + e — ^/a 2 + e 2 — 2ae + 4d/ has negative real part if and only if 



a 



+ e > ^a 2 + t 



that is, (S M + Me + A 7 ^)( Mm + a M A + A 



A I 



1 - 



> o 



2ae + Adf 
:)(rS M (l-e)) 



-) 



l-Rl > if fl| < 1 



where Rl. = V -LJ^ n f^tt- The term ~ gives the ratio of the infcc- 

r<5 (1 — e) 

tion to the total mortality rate of the erythrocytes while —. — - — ^ is the rate of release of 

merozoites against drug and immune action per dying merozoite . This means that R\. is the within-host 
reproductive potential of a merozoite released in presence of helminths larvae, against drug and immune 
action. Thus, the following conclusion is made: 

Lemma 2.2.3 The initial invasion state, the immune equilibrium where T = £:z r ± is stable if 



a rS — (1 



{K + Me + A/ 



£ ^)(M 



M 



< 1. 



(29) 



This biologically important criterion determines how malaria parasites or helminths larvae following a 
successful invasion, will eventually be controlled by the immune response. If the criterion is not satisfied, 
neither the malaria parasites nor the helminths larvae will be eliminated but either or both will persist 
at an equilibrium density within the host. Alternatively, if a host undergoes successful therapy of either 
infection or both, and the immune cells have settled at their equilibrium level, this criterion determines 
whether a new infection is able to overcome the stability of this equilibrium and re- infect the host. The 
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outcome of infection depends upon the balance between malaria parasite and helminths larvae and host 
attributes. As the death rate of the merozoites increases, the immune equilibrium becomes more and 
more stable and the host is able to eliminate the infection more easily. Similarly, the higher the equi- 
librium level of immunity p " b A ' t , the greater the malaria parasite invasion rate a M must be to achieve 
reinfection. It is interesting to note that in this system, the stability of the immune equilibrium does not 
depend upon the rates of T cell activation by the pathogens. However, the speed at which the system 
moves towards equilibrium will be influenced by the magnitude of these rates. 

For global stability of the initial state equilibrium with immunity, consider the following Lyapunov 
function, 

C = I e +M + H + T 
which is positive definite V I e , M, H, T as t — > oo, with orbital derivative given by 



£ = I' e + M' + H' + T' 

= a M U e M + a H U e H - [i e I e - S M I e - X : I e T + rS M (1 - e)I e - fJ, M M 
-a M U e M - X M MT + U- f M H H-5 H H-a H U e - X H HT + ^ M M 
+ ll I e + la H+pT - fi t T - bT 2 

= -Heh - S M I e + rS M (l - e)I e - [i M M + II - fi H H - 5 H H - a H I e H 
+ ( 7m - \ M T)M + ( 7h - \ H T)H + ll I e - \jI e T + P T- Mt T - bT 2 

< rS M (1 - e) + 7/ - (Me + K + X,T)I e + (p - H - bT)T 

= (Me + S M + X^ C^-^T) IK. + ^ 2 (^ - 1) 

- (Me + S M + X T T)(R 2 l)I e + bT 2 {^ i) whcrc R 2 = 

< if R 2 < 1 and since ^ < 1 



Ri is defined as the within-host rupture potential of an infected erythrocyte with immune action. Thus, 
the following result is obtained: 

Lemma 2.2.4 In presence of immune response, the initial invasion state is globally asymptotically stable 
in T if R 2 < 1. In this case, the immune population has been triggered but the malaria parasites and 
helminths larvae have been eliminated and the system has returned to an infection-free state with a degree 
of residual immunity. 

This further leads to the following conclusion: 

Lemma 2.2.5 It is further necessary that to avoid reinfection from malaria, p ~ b fJ ' t <C 1 since the higher 
the equilibrium level of immunity p ~,^' , the greater the malaria parasite invasion rate a M in presence of 
helminths larvae must be to achieve reinfection. 
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When Rf > 1, it implies that 



rS M (l-e)+ 7l ) 
(p, c +8 M +\ I T) ^ 



l-e> ^■(^e+^+A J T)-7 J 
e < 1 -77-(Me + ^M+A J T)+ 7i 



This expression gives the critical level of drug action below which no recovery is expected. Thus we have 
the following result: 

Theorem 2.2.1 The percentage threshold value of the drug action that is sufficient for recovery in pres- 
ence of immune response, ej is given by 

e i = l--±-(» e + 6 M +\ 1 T)+ % (30) 
Note that ej < e . Conversely, the following counter result is achieved: 

Theorem 2.2.2 The maximum value of T -cells required to neutralize the antigens to eliminate a co- 
infection of malaria and helminthiasis in presence of therapy, T max is given by 

T m ax = (rS M (1 - e) + 7, - /i e - 6 M ) (31) 

For endemicity, from (29), when T = ^p- we have cd - (a + XiT){b + X M T) > where a = S M + 
fi e , b = fi M + a M ^-, c = a M and d = rS M (1 - e). This gives (a + A/T)(6 + \ M T) - cd < implying 
that XjX^jT 2 + (b\ I + a\ M )T + (ab - cd) < and hence, f(T) = AT 2 + BT + C < 0. At steady 
state, the non-zero response of immune response satisfies f(T) = AT 2 + BT + C = 0, and this can be 
analyzed for the possibility of multiple immune responses. Note that the coefficient A = X T X M > and 
B = A J (/i M + oe M -^-) + X M (5 M + yu e ) > 0. Hence the following result is established: 

Theorem 2.2.3 There exists 

(i) precisely one positive immune response if C — [(S M + fi e )(fi M + a M -^-) — a M rd M -^-(l — e)] < 0, 

(ii) precisely one unique immune response if B 2 — AAC = [\j(p M + a M ^ ) + X M (8 M +^e)] 2 - i\ I X M [(8 M + 

(Hi) precisely two immune responses if B 2 — AAC > 0, 
(iv) The immune response is impaired otherwise. 

The following proposition states the local stability of the endemic equilibrium. 
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Proposition 2.2.1 The endemic equilibrium is locally asymptotically stable whenever it exists. 
Proof: 

At endemic equilibrium, from Jacobian 26, it is noted that since the diagonal entries of the Jacobian J e 
are negative and the eigenvalues of any square matrix are the same as those of its transpose, then using 
Gersgorin discs [26], it suffices to show that J e is stable if it is diagonally dominant in columns. Let 
fj, = max{gi, g 2 , g 3 ,g4, 35} where 

gx = -iJLe-a M M* -a H H* + a M M* +a H H* - a M M* -a H H* 

= -n e -a M M*-a H H* (32) 
< 



.92 = S M -n e - A/T* - a M M* + r5 M (1 - e) - a H H* + 7/ 

= —o. M M* - a H H* + ( Me + 5 M + XiT*)(Rf - 1), 

< 

33 = -Q-m Ue + a M U* - M M - a M U* - X M T* + 7m 
= -Mm - a M U* - \mT* + 7 M 

< since m m > 7m! (refer to Hcthcote, 2000). 

34 = -a H U* + a H U* -fi H -a H U* - X H T* + 7 H = (33) 
= -Mh - a fl (U* +1*)- ^ T* + 7 H 

< since fi H > j H 

35 = -A//* - A M M* - \ H H* + P -^ t -2bT 
= -XiI*-X M M*-X H H*-2bT + P -n t 
= -A 7 7* - A M M* - A H ff* - 2KT(1 - ^) 

< since T > from equation (25). 

^? = ^~+|~ +1j"t* ^ s * ne basic infection induced rupture rate with immune activation per dying erythro- 
cyte, then, Rf < 1 for stability. This implies that 

M = max{-/i c - a M M* - c^iT, -a M M* - a H H* + (/i e + 5 M + A/T*)(i? 4 2 - 1), 

-Mm ~ «m ~ A M T* + 7 M , -Mh - " H CC - A ff T * + 7ff i (34) 
-X T r e - X M M* - X H H* - 26T(1 - ^* ) 
< 0; hence diagonal dominance and conclusion of the proof. 

From the proposition, the following conclusion is made: 
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Lemma 2.2.6 The endemic equilibrium is locally asymptotically stable if Rf = — ^ +\jT* > ^ la ^ 
if the basic infection-induced rupture rate per erythrocyte with immune activation is less than the total 
mortality rate of the erythrocytes. 



3 Estimating the severity of the co-infection 

Let the basic acting factors of a co-infection of malaria and visceral helminthiasis be concentration 
of pathogenic multiplying antigens, H(t), and M(t), defined in the previous section, concentration of 
antibodies, A(t) = A M + A H , where A M and A H are the antibodies against the malaria parasites and 
helminths larvae respectively; concentration of plasma cells C(t) and the relative characteristic of the 
affected organ, m(t), in this case, red blood cell depiction. Consider the following model:- 



^ = (f3 M - lM A M )M 

f: = W h - 1h a h )h 

= a M A M M + a H A H H-fi c (C-C°) 

= p M C-fi A A M -r] M ^ M MA M 

= p H C- f i A A H - VHlH HA H 

% = u M M + a H H-n M m 



where (3 M is the average rate of malaria parasite multiplication, and j3 H is the average rate of helminths 
larvae invasion of red blood cells, j h AhH and j m AmM is the number of antigens neutralized by the 
respective antibodies during the time interval At. 7 M and -f H are therefore parameters connected with 
the probability for an antigen to be neutralized after encountering the antibodies. The third equation 
describes the growth of plasma cells. This is the population of plasma cells, producers of antibodies; in 
this model, the antibodies are the substrates capable of binding the malaria parasites, and the helminths 
larvae, therefore, the number of lymphocytes stimulated in this way is proportional to AmM and AhH. 
This is the relationship describing the increase of plasma cells over a normal level C°, that is the constant 
level of plasma cells in a healthy human. a M , a H are defined as the coefficients allowing for the probability 
of antigen-antibody collision respectively, the stimulation of cascade reaction, and the number of newly 
generated cells. We define as the time it takes the plasma cells to age, implying that p c (C — C°)At 
is the decrease in the number of plasma cells due to aging [28] . 



Let p M C and p H C be the generation of antibodies fighting malaria parasites and helminths larvae by 
plasma cells, giving p = p M + p H as the total rate of antibodies production by the plasma cells. The total 
number of eliminated antigens that were neutralized by antibodies is equal to ^ m AmM + j h AhH, so if 
neutralization of one antigen requires r\ M +rj H antibodies, then this implies that t] m j m AmM + iti h ^i h AhH 
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describes the total decrease in the number of antibodies due to binding with malaria parasites and 
helminths larvae. The total drop in antibodies population due to aging, where \x A is the per capita rate 
of decay of antibodies, is given by \i a Am + [i a Ah- 

Next we consider depletion of the red blood cells during the co-infection of malaria and helminths, 
due to the decrease in the activity of organs providing the delivery of immunologic materials, such as 
leukocytes, lymphocytes, and antibodies, that are necessary for the struggle with multiplying antigens. 
Assuming that the productivity of the red blood cells depend on the degree of invasion by malaria parasites 
and helminths larvae. Let R be the relative characteristic of a mass of healthy red blood cells and R d 
the corresponding characteristic of the healthy red blood cells in an invaded red blood cell mass. Define 
m as 

1 - TP that 

is, m is relative characteristic of damage to the red blood cell mass. For non invaded 
red blood cells, m is equal to zero naturally, and for completely impregnated red blood cells, m equals 
1. cr M M + (7 H H signifies the total degree of depletion of red blood cells. It is assumed that the increase 
in relative depletion of red blood cells is proportional to the number of antigens M + H , described by 
<J M M + <7 H H, where <r M , a H are constants. A decrease in damage is caused by the recovery activity of 
an organism, that is, fiM is the restoration of the red blood cells. 

3.1 Steady states during the severity of the co-infection 

For stationary solutions, we equate the derivatives to zero and solve with initial conditions. Note that 
at t = 0, C = C°, A* M = Am, A* h = Ah are the values for a non invaded red blood cell mass for 
M = = H . This implies that the trivial solution, describing the healthy state of red blood cells is 

M = 0, H = 0, C = C°, A* M = ^-,A* H = m = (36) 

Ma Ma 

Thus we have the following conclusion: 

Lemma 3.1.1 The concentration of antigen population and damaged mass of red blood cells at disease- 
free state are equal to zero, and the quantities of plasma cells, C and antibodies A M and A H correspond 
to the values of an immunological status of a healthy individual. 

Consider the system of points (M, H, C, Am , Ah , m) in the neighborhood of the equilibrium point 
(M*,H*,C*,A* M ,A* H ,m). Let M(t) = u{t), H(t)=v(t), C(t) = C*+w(t),A M (t) = A* M +x(t),A H (t) = 
A* H +y(t), m(t) = z(t) where u, v, w, x ,y, z are small deviations of functions from this equilibrium state 
(36). Substituting these expressions into the system of equations (35) and assuming that the deviations 
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are small, neglecting small variables of the second order of smallness we obtain: 



f = (/3 M - lA* M )u 
f t =(P H -lA H )v 

W = a M A *M U + a H A ll V -Vc W 

= Pm w ~ Pm1m A *M u - Pa x 
% = Ph w - Ph1h A *H v - PaV 



<J„V 



Pmz 



(37) 



The corresponding Jacobian 
/ 



dJ 
~dt 



Pm 1m A *M 




















Ph ~1h A H 














a M A M 


a H A H 


-V-c 











Vm 7m a *m 





Pm 


-Pa 











-Vh1h A *H 


Ph 





-Pa 





CT M 















z) T (t), with 


distinct eigenvalues 


— PM, 


~Pa,~ 





(38) 



Pm J 



result: 



Lemma 3.1.2 All small deviations from stationary solution (36) for j3 M < j M A* M and j3 H < j H A* H 
tend to zero with time, that is, the total number of malaria parasites multiplication, plus the total number 
of helminths larvae that invade the red blood cells at a given time must be less than the total number of 
malaria parasites and helminths larvae neutralized by the antibodies during an episode when both malaria 
parasites and helminths larvae are present in the red blood cells. This implies asymptotic stability of the 
invasion-free state. 

This leads to the following conclusion: 

Lemma 3.1.3 From (35), at steady state, the small infection dose that does not lead to the loss of 
stability is estimated by 



° fl---v..Ti ' o H 7 H n H 



J M 'M 'M 



(39) 



where M* and H* are the immunological barrier values for malaria and helminthiasis respectively. 



The immune barrier is exceeded when an infection dose M Q and H Q satisfies conditions M Q > M* and 
H Q > H* , and is not exceeded otherwise. Biologically, it implies that, if in the case of red blood cell 
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invasion by a small dose of malaria parasites and helminths larvae, the immunological barrier cannot be 
exceeded, then regardless of the infectious dose, the invasion does not develop, that is, the number of 
malaria parasites and helminths larvae in the organism decreases with time tending to zero, and the red 

blood cells are restored. In addition, the elevation of C°, that is of the level of immunocompetent cells in 

p c° p c° 

a healthy red blood cells mass increases the immunological barrier, since A^ — ~~ ~ — ~ an d A^ — — — — — 
and therefore, it is an effective method of prophylaxis, and possibly, of disease treatment. This leads to 
the following lemma. 

Lemma 3.1.4 The immunocompetent cells threshold values against malaria and helminths are given by 
0<M o < 1 ^ c °m- 1 =M * o < H < ~IhPhC° h -i , , . 

with C° M — - 1 p and C° = - 1 p as the threshold values of the immunocompetent cells respectively, 
below which, no immune response is expected. 



3.2 Results of the severity of the co-infection 

Chronic form of either disease is caused by an insufficiently effective reaction of the immune system 
that happens when a M or a H is small, that is if the immune system's reaction is weak. In this case, 
M — > M max and H — > H max ; and on the other hand, A* M and A* H drop respectively. The outcome of the 
co- infection depends on the derivatives and ^ , if they can become negative and how long they do 
not change signs. It is noted that ^ < if M(t) > and A M (t) > ^ and ^ < if H(t) > and 
A H (t) > ■ Thus we have the following conclusion: 

Lemma 3.2.1 If we assume that infection of either or both malaria and visceral helminthiasis has hap- 
pened, then the A M (t) > and/or A H (t) > is the necessary and sufficient condition for and 
to be negative. 

If the immunological barrier is not to be exceeded, (M D < M) and (H a < H), then only CASE I of the 
following four cases is possible:- 

CASE I ^M- < 0, =§r < in a large interval of time. In this case, there will be a subclinical form of the 
co-infection. This happens since f3 M < "1 M A* M and j3 H < j H A* H in which case either there is an effec- 
tive (normal) immune response, when a M (p M + p H ) > PcVmJm an d a H (PH + Ph) > PcVhIh? an d 
weak response, (immunodeficiency) when a M (p M + p H ) < Pc'HmIm an d a H (PH + Ph) < PcVhIh- 
That is, when the infection dose is small, lower than the immunological barrier values, the re- 
moval of the malaria parasite or helminths larvae from the red blood circulation depends neither 
on the doses of infections nor on the strength of immune response. The removal of either antigen 
is provided by the antibody levels A* M or A* H against each disease that are present in the body at 
the time of infection. When the doses of malaria parasites and/or visceral helminthiasis increase 
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to considerable extents as compared to the immune barrier, the strength of the immune response 
begins to play an important role. 

In case of fi M > j m A* m , (3 H > j H A* H , ^M. > 0, 4M. > for an interval (ti,t 2 ), and assuming that 
M — > M max and H — > H max at a point ti in the interval, and then decrease to M ^ M min and 
H — > H m i n , we have two cases: 
EITHER 

CASE II ^ < 0, < within a sufficiently large time interval (£1^2) where t 2 is the time when 
M — > M mi „ and ff — > H min . Hence, acu£e /orm 0/ i/ie co-infection is obtained in which case, f3 M > 
1m 7 ! mi Ph > IhVh an d therefore the immunological barrier against malaria parasites or helminths 
larvae does not exist. Therefore, there would be a rapid increase of malaria parasites and helminths 
larvae in the red blood cells up to the values exceeding the infection doses, and then rapidly get 
eliminated due to a strong and effective immune response that leads to the production of antibodies 
against the malaria parasites and helminths larvae in quantities sufficient for the elimination of both 
infections. Note that the higher (3 M and f3 H are, the higher A* M and A* H respectively; hence the 
faster the maximum value of malaria and helminths fighting antigens quantities are attained and 
the faster the process will stop. This gives the following result: 

Lemma 3.2.2 With a high rate of malaria parasite multiplication, or a high rate of red blood 
invasion by the helminths larvae, or more generally, a high dose of infection by either disease, with 
other parameters constant, the amount of malaria parasites and helminths larvae that stimulate the 
immune system effectively is attained more quickly. 

This implies that M max or H max does not depend on the dose of infection but is determined by 
the immune status in relation to the malaria parasites or helminths larvae; that is, by the set of 
model parameters. The dose of the infection affects the time when M max or H max is reached. In 
this form of infection, there is red blood depletion due to malaria, a M and due to helminths, a H . 
OR 

CASE III ^ < 0, ^ < within a sufficiently small time interval (t\, t 2 ), then the chronic form of the 
co-infection is expected. 

If however a point t\ does not exist, that is if, 

CASE IV ^ > 0, > within an infinitely large interval of time, then a lethal outcome of the 
co-infection is obtained. This happens when the formation of plasma cells is delayed and in turn a 
delay in the production of antibodies specific to malaria parasites or helminths larvae. This leads 
to the following result: 
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Lemma 3.2.3 In order to prevent a hyper-toxic form of the co-infection of malaria and helminthi- 
asis, it is necessary to delay the replication rate of malaria parasites (3 M and the invasion rate of 
helminths larvae, 13 H . 



3.3 Threshold values for the severity of the co-infection 

If both or one of the antigens is eliminated, then from (35) at steady state, the following equilibria exist: 



E 
Eh 



E 



M 



= (0,0, c,£ 
= 



,c p h c 



Pa ' Pa 

P A P C 



, 0) the disease-free, 

a H f } H^A 



(a 



hPh-Pc^h^h) 
PaPc 



mPm Pc^m^m") 
P A Pr 



0. 



a M p M PJ 



H M ' M 



0» 



(41) 



where endemicity of either equilibrium depends on the respective threshold values. 

For competitive dominance of helminths larvae by the malaria parasites, the following equilibrium point 
will be attained : 



E 



M H 



Pa 



P M Pjj °mPa 



(Pm+Ph)+Pc^m (Pm 



(42) 



with R? M > 1 > R 2 Ho where R eM is the basic reproductive potential for a successful malaria dominance 
in presence of helminths and 

R 2 _ a M (P M +Ph) n2 _ Pc^hIh (43) 

Thus, we obtain the following result: 

Lemma 3.3.1 For competitive exclusion of visceral helminthiasis by malaria, R > 1 > R H . That is, 
the total rate if malaria parasites- antibody collision must be greater than the loss due to neutralizing by 
the antibodies; in this case, the malaria parasites have become resilient in their quest for survival. 

If helminths dominates malaria, then the endemic equilibrium will be given by 



W TT _ PaPc Pa Pa i a H (f> M +Ph )+Pc Vh (£h ZlM )] P M 

Hm <*m(Pm+PhK1-K 2 oM )' V H 7 H (Rl H -l)> PhPc^h1h{Kh-^ ' 7« ' (44) 

Ph. a H PA 
1 H ' CMiJ H 7 H (^ H -l) 
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which exists whenever R 2 H > 1 > R? M where R eI1 is the basic within- host invasion potential of helminths 
larvae in presence of malaria parasites and 



?2 _ "h(Pm+Ph) 

This gives the following similar result: 



£2 _ '±h}Tm_J±h. ) ft2 _ Pc^m^m 
h M C 1 H 7 H »« "m^m+Ch) 



(45) 



Lemma 3.3.2 If visceral helminthiasis dominates malaria, R 2 > 1 > R 2 . Similarly, the helminths 

•> ? eH oM 07 

larvae are resilient but the malaria parasites are being neutralized by the antibodies, and so R 2 oM < 1 ■ 
However, if both antigens compete for coexistence, then one of the following endemic equilibria is expected: 
A. Case of weak sub thresholds 

CASE I RoM < 1, ReH > 1, ReMH > 1 



E H M = 



( P A Kf [Pm +Ph )+Pq 1m ZZm )] 4. Pa l a H (Pm +P H )+Pc V h [£h Zln )} 

V PM a M(PM+P H )( 1 - R2 oM ^ 



PHPcriH-tH&in)- 1 ) 
f> M Ph \ a M Pa P% Tm a M (Pm +P H ^'Km^h Pa Pq 11 h ■ h (_Kh 



Tm ' Th ' 



PM ( Q M a H (PM +P H ) 2 +P 2 c Vm V H 7„ Tff ) ( R l 



CASE II R e M > 1, RoH < 1, ReMH > 1 



(46) 



E m h = 



%7 M (^ H -l)' Qfl (p M +P H )(l-^ H )' 



j mPc 



+ 



,P A P_c2mJjA R , 



Pm. Eh. 

7 M ' 7„ ' PM(a M a H (p M +p H y+^ cVMVH t M j H )(Rl MH -l) 



Ph°'h(Pm+Ph)( 1 ~ R2 oH ) 

-i)+<y H n A n 2 c n H i H Ci--R 2 oH )] 



B. Case of strong sub thresholds R m < 1, RoH < 1, RoMH < 1 



(47) 



HM 



PaP 



_ PaPc 

r ("m +Ph j(l-B^ M ) ' ~n (Pm +Ph)^- r I h 5 ' 
f Kf (Pm +Pff )+^c (^m ~Tm )1 i ["g (Pm +Pff )+Pc r h i %j Z2m )] ^ 

/3 M /3 H [ g M A 1 A *7m T m "m (Pm +Ph ^'Km }+° 'h Pa P% r >n'<H a H (Pm +Ph K___Kh}\ 



Tm ' T H ' 



P™P C (PM+PH)( a M r >H'rH+ a H r >M'rM)( 1 - R2 nMH ) 



(48) 



where R 



Pc (Pm +Pg ) ["m^h Tfj + a g Tm t ?m 1 J and 
[«m "h (Pm +Pff ) 2 +^c Iff Tm Th ] 



[ q m a ff (Pm +Ph ) 2 +Pc ''m »?h T m T„ ] 



oM H - M C (P M +Pff ) K-f J '« Tff +"„ 7 M f M I 



< 1 is the within-host basic reproductive potential for a success- 
ful co-infection of malaria and helminths. This gives the following result: 
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Lemma 3.3.3 The malaria parasites and helminths larva will co-exist, for strong sub thresholds and 
one of them will dominate the other for weak sub thresholds. This implies that there may be competitor- 
oscillatory coexistence until the system bi-stabilizes to endemic equilibria, depending on the immune action 
against either or both antigens. 

Further from the above lemma, the following conclusion is made: 

Lemma 3.3.4 For endemicity of the within-host co-infection of malaria and visceral helminthiasis, one 
of the reproduction numbers must be greater than 1 and coexistence equilibria exists; no coexistence exists 
when both reproductive potentials are greater than 1. 

4 Discussion 

In the model for the co-infection of malaria and visceral helminthiasis, we have assumed that an individ- 
ual gets malaria infection during invasion by helminth larvae. In absence of immune response, our results 
show that there is no disease-free state, but an endemic state is attained via damped oscillations towards 
a fixed point (see Figure 1 A,C and B). Our results showed that antigens invade the blood system if 
the rate of red cell rupture per invading merozoite is greater than one as seen in Lemma 2.1.1. If more 
merozoite are released, there exists malaria-only endemic point. However, both antigens co-exist if the 
mean infection burden is greater than one. In this case, there is a threshold value for drug action below 
which no recovery of host is expected, Lemma 2.1.4. In presence of immune response, three equilibrium 
states exist. The initial invasion state, and secondly, the unstable state when the immune population 
has been triggered and the antigens have been eliminated. The third state represents the endemic state 
which is stable if the infection-induced rupture rate per erythrocyte with immune activation is less than 
the total mortality of the erythrocyte (see Lemma 2.2.1 and Figure 2). A model for severity of the 
co-infection shows in Figure 2B that an immune response will be delayed until immunological barrier 
values for malaria and helminth are exceeded. Our results show further that once the immune response 
is triggered, the system in inherently unstable and will not return to zero, Lemma 2.2.2. 

In this paper, a co-infection of within-host malaria and visceral helmnithiasis is studied. The interest 
was to study immune response and impairment when bot antigens invade the red blood cells. The process 
of immune response is more complicated than modelled here, but the model has important consequences 
for the success of disease control programs. This model is a conceptually simple and a logical extension 
of the classic individual species models that have formed the foundation of epidemiological theory over 
the last 30 years. As such, it should provide a valuable framework for understanding the dynamics 
of parasite co-infection with immune response and evaluating the wide-ranging implications of disease 
control strategies in the future. Although with constraints due to limited data, the study shows that 
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the epidemiology and population dynamics of within-host-microparasite-macroparasite communities are 
highly complex, resulting in different predictions from those made by the corresponding single parasite 
species models. If better estimates of parameters become available, more detailed studies can be done by 
taking into account some of the aspects that have been overlooked. 
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